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CN ' We extend the quantal hypernetted-chain (QHNC) method, which has been 

a^. 

Q^ , proved to yield accurate results for liquid metals, to treat a partially ionized 



plasma. In a plasma, the electrons change from a quantum to a classical fluid 
gradually with increasing temperature; the QHNC method applied to the elec- 
tron gas is in fact able to provide the electron-electron correlation at arbitrary 



^ ' temperature. As an illustrating example of this approach, we investigate how 



liquid rubidium becomes a plasma by increasing the temperature from to 
30 eV at a fixed normal ion-density 1.03x 10^^/cm^. The electron-ion radial 
distribution function (RDF) in liquid Rb has distinct inner-core and outer- 
core parts. Even at a temperature of 1 eV, this clear distinction remains as 



a characteristic of a liquid metal. At a temperature of 3 eV, this distinction 
disappears, and rubidium becomes a plasma with the ionization 1.21. The 
temperature variations of bound levels in each ion and the average ionization 
are calculated in Rb plasmas at the same time. Using the density-functional 
theory, we also derive the Saha equation applicable even to a high-density 
plasma at low temperatures. The QHNC method provides a procedure to 
solve this Saha equation with ease by using a recursive formula; the charge 
population of differently ionized species are obtained in Rb plasmas at several 
temperatures. In this way, it is shown that, with the atomic number as the 
only input, the QHNC method produces the average ionization, the electron- 
ion and ion-ion RDF's, and the charge population which are consistent with 
the atomic structure of each ion for a partially ionized plasma. 
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I. INTRODUCTION 

In order to calculate thermodynamic functions, transport coefficients and optical prop- 
erties in a partially ionized plasma, it is a fundamental problem to determine the average 
ionization Zj, the equilibrium correlations among ions and electrons, the atomic structure 
(bound-levels in the ions) and the charge population (ionization balance) of differently ion- 
ized species, in a self-consistent way with each other in these quantities. However, at the 
present stage there is no theory which can produce these quantities in a partially ionized 
plasma in a unified manner. It is the purpose of the present work to show that the quantal 
hypernetted-chain (QHNC) equation [||] developed for a liquid metal can be extended to 
calculate these quantities of a partially ionized plasma in a unified manner. 

Up to the present, in the calculation of thermodynamic functions or optical properties in 
a partially ionized plasma, the ion-sphere (IS) model [@-|^ is used as the standard method. 
Although there are many kinds of variations in the IS model, the essential point of this model 
is that the ion-ion correlation in a plasma is approximated by the step function 9{r—a) with 
the Wigner-Seitz radius a; an atom is considered to be either confined within the ion-sphere 
or immersed in the infinite jellium. It should be remarked that this model is only applicable 
for high-density and low-temperature systems, that is, limited to a narrow range of densities 
and temperatures for the plasma state. 

On the basis of the density functional (DF) theory, Dharma-wardana and Perrot P] 
derived a set of integral equations for static correlation functions in a strongly coupled 
plasma; it was shown P| that this method breaks down when it treats a plasma with a 
significant number of bound electrons as an ion due to its improper treatment of the electron- 
ion correlation. To judge whether a theory treating a partially ionized plasma is proper or 
not, we have the criterion as to what extent it can reproduce the observed structure factors 
of liquid metals for which many reliable experimental data exist. This is a liquid metal can 
be taken as a special type of partially ionized plasma. In this context, we have proposed a set 
of integral equations for radial distribution functions (RDF) in a liquid metal as a nucleus- 



electron mixture |]l| on the basis of the DF theory in the quantal hypernetted-chain (QHNC) 
approximation. At this point, it should be kept in mind that the QHNC theory is limited to 
treating only a "simple" liquid metal, where the bound states are clearly distinguished from 
the continuum state. Already, we have applied the QHNC method to several simple liquid 
metals |Ty-|I^ , and obtained their structure factors in excellent agreement with experiments. 
In these calculations, we have demonstrated that the QHNC method can determine the 
"outer structure" (the ion-ion and electron-ion RDF's and the ionic charge Zi) in a consistent 
way with the "inner structure" (the atomic structure of the ions) using the atomic number Za 
of the system as the only input data. Therefore, the QHNC method is suited for treating 
a plasma, where the ion-ion and electron-ion interactions may vary over a wide range in 
conjunction with the internal structure of each ion according to change of state condition. 



In a similar spirit, Perrot ]T5| has proposed the neutral-pseudoatom (NPA) method based 
on the DF theory to calculate the effective ion-ion potential for a partially ionized plasma. 
Gonzalez et al. have successfully applied this theory to alkali liquids [ITB] and alkaline-earth 



liquids |T^. The NPA method can be derived from the QHNC theory with additional use 
of the IS approximation |^ in the determination of the pseudopotential to construct an 
effective ion-ion interaction. This model is called in our approach the jellium- vacancy model 
I nil and is used to obtain an initial guess for the effective ion-ion interaction in the iteration 



to solve the QHNC equation for a liquid metal. At this point, it should be recognized that 
the NPA method is not appropriate for treating a high-temperature plasma with the weak 
ion-ion correlation, since this method is based on the IS model; this fact will be discussed 
in the present work. 

From successful applications of the QHNC method to liquid metals as mentioned before, 
we can expect this method constitutes an adequate representation of a partially ionized 
plasma. In the present work, we extend the QHNC method applicable to a plasma state, 
where the electrons change from a quantum fluid to a classical fluid gradually with increasing 
temperature; this is in contrast with a liquid metal where the electrons can be assumed to 
be perfectly in the Fermi degenerate state at zero temperature because of the density being 



high. When the electrons in a plasma begin to behave as classical particles, it is difficult to 
calculate the free-electron density distribution under the external potential by solving the 
wave equation. Therefore, Furukawa |T^ used the Thomas-Fermi (TF) approximation to 
evaluate the electron-ion RDF in the QHNC equation for a plasma. Also, Xu and Hansen 
TP| applied the TF version of the QHNC method to a hydrogen plasma with a gradient 



correction to the TF kinetic energy of the electrons. In the full use of the wave equation 
to calculate the free-electron density distribution, we demonstrate that the QHNC method 
can treat a partially ionized plasma by taking liquid rubidium as an illustrating example; 
this exhibits what changes are found when a liquid-metal turns to a plasma state with 
increasing temperature at a fixed ion-density. In a liquid metal, there is clear distinction 
between the inner-core and outer-core structures in the electron-ion RDF, which allows one 
to construct an electron-ion pseudopotential in a liquid metal. When a liquid metal changes 
into a plasma state, this distinction disappears, and we can not set up a peudopotential in 
the same manner as is performed in the usual liquid-metal theory. 

In principle, the DF theory generates the exact density distributions of electrons and 
ions in a plasma. Note that it can yield only the average-ion structure in a partially ionized 
plasma. In a real plasma, there are many differently ionized ions around the average ion. The 
fundamental theory of ionization in a plasma is provided by the Saha equation. However, 
the usual Saha equation can be applied only to a low density equilibrium plasma, where the 
interactions among particles are negligible. Although there are many modifications of the 
Saha equation for applicability to dense plasmas by introducing the continuum lowering, any 
modified theory can not treat a dense plasma at low temperatures. In the present work, on 
the basis of the DF theory we derive the Saha equation, which provides the charge population 
of differently ionized ions in a dense plasma in the region from low to high temperatures; this 
charge population yields an average ionization consistent with an average ion in the plasma 
determined the DF theory. The bound levels and the chemical potential contained in the 
Saha equation are supplied by the QHNC equation for the ion-ion and electron-ion RDF's 
in the plasma. In this way, the QHNC method is shown to generate the electron-ion and 



ion-ion RDF's, the average ionization Zj, and the charge population of differently ionized 
species to be consistent with the atomic structure (bound levels) of each ion in a unified 
fashion, as will be shown by the example of a rubidium plasma. 

The paper is organized as follows. In §. |I|, we present a summary of the QHNC method 
along with the one-component QHNC equation for an electron gas in the jellium model, 
making extensions to treat a plasma. As an illustrating example, the application of this 



formulation to a Rb plasma is shown in § |T|, where the numerical technique to solve the 



QHNC equation is explained. In § |^ we set up the Saha equation on the basis of the 
DF theory, and the charge populations are calculated for Rb plasmas at a fixed liquid- 
metal density for several temperatures. The last section is devoted to discussion, where 
the limitation of the IS model is also examined, and prospects of applications based on the 
QHNC method are mentioned, such as calculations of the atomic structure and transport 
coefficients in a plasma, and the determination of the effective interactions to be used for 
the molecular-dynamics simulation of a plasma as a classical electron-ion mixture. 

II. LIQUID METAL AND PLASMA AS NUCLEUS-ELECTRON MIXTURE 

We can set up the following three models to treat liquid metals and plasmas. In the 
simplest model, a liquid metal can be treated as a neutral liquid, when its interatomic 
potential is given. However, this interatomic potential must be introduced from the outside 
of the model; this is impossible when treating a plasma. In the second model, a liquid metal 
is considered as a binary mixture of ions and electrons; in this model, the ionic charge Z\ and 



the electron-ion interaction are unknown except for a perfectly ionized plasma pO|J^, even if 
the ion-ion interaction is taken as a pure Coulombic. In the case of a plasma, it is rather a 
fundamental problem to determine the ionization Zi. Most fundamental model is to consider 
a liquid metal as composed of nuclei and electrons. In this model, all input data are known 
beforehand if provided the atomic number Z^ of a liquid metal; this first-principles approach 
enable us to treat a plasma in a wide range of temperatures and densities. Therefore, let 



us think of a plasma as a nucleus-electron mixture |l| consisting of A^i nuclei and Zp^Ni 
electrons. Here, we single out one nucleus and fix it at the origin; then, the fixed nucleus at 
the origin in the mixture causes the external potentials for electrons and ions, and induces 
an inhomogeneous system. This inhomogeneous system can be equivalently translated into 
a simpler system: a fixed nucleus with the atomic number Z^ is surrounded by electrons 
and ions, of which structure pb(^) is undetermined at first and should be determined self- 
consistently at final. In this simplified model, the fixed nucleus at the origin is surrounded 
by (iVi— 1) ions interacting via a potential fii(r) and by Zi{Ni—1)+Za electrons. Under this 
circumstance, the DF theory provides effective external potentials v1^{r) {i=e or I), which 
yield the exact electron- and ion-density distributions around the fixed nucleus on the basis 
of the following reference system [^: a mixture composed of (A^i — 1) noninteracting ions 
and Zi{Ni—l)+Zji^ noninteracting electrons. Each ion in the reference system is assumed to 
have Zb bound-electrons with a distribution ph{r). The electron-density distribution around 
the central nucleus is determined by solving the wave equation for noninteracting electrons 
in the form: 

ne(r|N) = n°(r|t;,^^) = n,^(r|N) + n^(r|N) , (2.1) 

under the effective external potential, 

vtUr) = -Z^e^r + -^^ - /.^ , (2.2) 

which is constructed by the DF theory in terms of the interaction part of the intrinsic 
free-energy jFj^t defined from the reference system and the interaction part of the chemical 
potential /x™*. Here, n^{r\U) is determined by solving the wave equation for an electron 
under the external potential U{r) in the form: 

nl{r\U) = J:Jie.)Mr)\\ (2.3) 

where 4'i{r) is the wave function of the state e^ (bound or free) and /(e), the Fermi dis- 
tribution function. Now, the density distribution ng(r|t>|^) can be divided into two parts: 



the bound- and free-electron density distributions, rig and nl, respectively, according to the 
states ej being negative or positive. 

The bound electrons n]^{r) =n^{r\N) around the nucleus at the origin should be taken 
to constitute the ion at the origin. Therefore, unknown ion structure phir) in the premise 
is determined by the condition that the central ion with n^{r) formed at the origin must be 
the same to any surrounding ion with the assumed structure Ph{r) around it, that is, 

Pb(r)=n^(r|N). (2.4) 

This leads to a self-consistent condition to determine the ion structure Phi'f')] the central ion 
structure n^{r) is iteratively used as the input ph{r) giving the new reference system, which 
determines the next effective potential (|2.2|) to evaluate a new n^{r). From this relation, the 
bound-electron number Z^ of surrounding ions can be evaluated from the bound-electron 
density distribution ng(r|N) by Zb = Jl^ n^{r)dr, that is, 

^^"Sexp[/3(e.-/.o)] + i' (2-5) 

for the ion with M bound states with the degeneracy Qi, and the ionic charge is obtained by 
Zi = Za — Zb- In addition, the chemical potential /Xg involved in Eq. ( |2.5| ) is determined by 
noting the fact: 

lim r.:(r|I) = / ,,p[^(^./2m - ^g)] + 1 (2!^)3 = "" = ^^^^ ' ^^.G) 

which can be rewritten in the form: 

7 _ V' 9i W 2 dp 

^~ti exp[/5(6, - /iO)] + 1 + 47 exp[/3(pV2m - p^J] + 1 {27rhr ' ^^'^^ 

Here, Uq denotes the number density of electrons or ions (z =e or 1). Also, the bare ion- 



electron interaction is obtained from Eq. ( |2.2| ) with use of some approximations ||2^ : 

v,,{r) ^ -^ + Jv:Mr - r'l) n^(r')rfr' + /ixcK^r) + nl) - /ixcK), (2-8) 

where p^ci^o) ^^ ^^^ exchange-correlation potential in the local-density approximation. 



With use of the average ionization Z\ and the bare ion-electron interaction fei(^), a 
plasma can be now modeled as a mixture of electrons and ions interacting through pair 
potentials Vijir) [z, j = e or I]. Applying the DF theory to this electron-ion mixture model, 
the ion-ion and electron- ion RDF's gair) are exactly expressed in terms of direct correlation 
functions (DCF) Cjj(r) and bridge functions Buif) as follows [Q: 

<?ei(r)=nO/(r|[/f)K, (2.9) 

(7n(r)=exp[-/3f/f(r)] , (2.10) 

with 

Uf{r) = v,^{r) - ra{r)/P - Ba{r)/P , (2.11) 

r,i(r) ^Y.1 Cu{\r - AW,[9ii{r') - l]dv' . (2.12) 

I ■' 

Here, n^J {r\Ulf) indicates the free-electron part of the density distribution. We can see 
a similarity of these expressions to those of a classical binary mixture. The difference is 
found only in that the electron-ion RDF is determined by a wave equation in stead of the 
Boltzmann factor exp[— /5t/|^(r)]. 

Moreover, it is shown from these expressions that the electron-ion mixture can be de- 
scribed as a one-component fluid interacting only via pairwise interaction feflf(?"), if the 
bridge function -Bii(r) is taken to be the one-component bridge function. Namely, the inte- 
gral equations ( p.9|) -( pTT^ ) for gii{r) can be transformed into a set of integral equations for 



the one-component model of plasmas. One is an usual integral equation for the DCF C{r) 
of a one-component fluid: 

C{t) = exp[-/3t;eff(r) + 7(r) + Bu{r)] - 1 - 7(r) (2.13) 

with an interaction Vesij), and the other is an equation for Vcs_{r), that is expressed in the 
form of an integral equation for the electron-ion DCF Cci(r): 

5Cei(r) = nV{r\v,i-r,i/(3-Bj(3)lnl - 1 - Br,^{r) , (2.14) 



since the effective interionic interaction Vcs{r) is given by 



l-n^oCociQ)x'Q 



f3vMQ) = PvuiQ) - J ':J.\^^^o ■ (2-15) 



Here, 7(r) = / C(|r — r'\)nQ[gu{r') — l]dr', and B denotes an operator defined by 

^Q[5"/(r)] ^ (xJ)"^Q[/(r)] = (xy"/exp[^Q ■ r]f{r)dr , (2.16) 

for an arbitrary real number a, and represents a quantum-effect of tlie electrons tlirougli tlie 
density response function Xq of tlie noninteracting electron gas. We can obtain a set of closed 
integral equations (referred to as the QHNC equation) from Eqs. ( p.l3|) -( p.l5|) by introducing 
the following approximations [|l|. (1) B^i — (the HNC approximation). (2) The bridge 
function Bu of the ion-electron mixture is approximated by that of one-component hard- 
sphere fluid (modified HNC approximation p3|). (3) An approximate Vei{r) is taken to 
be Eq. ( p.8| ), where we adopt the Gunnarsson-Lundqvist formula ||2^ for the exchange- 
correlation potential /^xc(^o)- (4) vii{r) is taken as a pure Coulombic Zle^ /r. (5) For liquid 
metals where the electrons can be assumed to be at zero temperature, the electron-electron 
DCF is approximated by Cee(Q) — — /9fce(<5)[l — G^'^^^iQ)] in terms of the the local-field 
correction (LFC) G^'^^^{Q) of the jellium model; we have used the LFC proposed by Geldart 
and Vosko pSf in our many applications to liquid metals. 

Under these approximations, a set of integral equations can be solved to determine 
the electron-ion and ion-ion correlations in a liquid metal together with the ionization and 
electron bound states. Figure |l] is an applied example of this procedure to liquid metals: 
the structure factors of liquid rubidium at temperature 313 K and density 1.03 x 10^^/cm^. 
The full curve is the QHNC result, which exhibits an excellent agreement with experiments 
denoted by open circles (neutron-scattering ||2^) and full circles (X-ray |^). From this 
example, we can expect the QHNC method to yield good results for partially ionized plasmas. 

However, when we apply the QHNC method to a plasma, there occur the following 
problems; the LFC involved in the electron-electron DCF must be evaluated at arbitrary 
temperature in dealing with a plasma. Although the LFC at the absolute zero temperature 
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has been calculated by many investigators and applied to many kinds of liquid metals, 
there is no standard way to calculate the LFC at finite temperature. For this purpose, we 



adopt the one-component QHNC equation |28] for an electron gas in the uniform positive 



background to obtain the electron-electron DCF at arbitrary temperature, which is written 
for an integral equation for the electron-density distribution ne(r|e) around the fixed electron 
in an electron gas, 

ne(r|e) = n°(r|[/eff) = S,/(e,)|V^,(r)|2 (2.17) 

with 

(3U,s{r) = (3vUr)- J C,Mr-rWe{r'\e)-nl]dr' . (2.18) 

Here, the DCF for a one-component system is defined by 

n^CUQ) = 1/X°Q - 1/Xq = -(3vUQ)[l - G{Q)] . (2.19) 

The Fourier transform of the density distribution yields the following bootstrap relation to 
determine the DCF Cee{Q) with combined use of Eqs. ( ^.IT] ) and ( |2.18| ): 



J^Q[n,{r\e) - nl] = Xq/x'q " 1 = l/[l-n^Cee(g)x^] - 1, (2.20) 



which is derived from a certain ansatz p8[. At this point, it should be noted that the 
QHNC equation (|2.17|) reduces to the well known HNC equation for a classical electron 
gas in the high temperature limit, because of the relations: Xq = 1; Xq = Sce{Q) and 
ng(r|?7eff) = nQexp[— /9f/cfj(r)] in the classical limit. As a result, the approximation (5) 
can be replaced by the one-component QHNC equation (|2.17|) for the electron gas; with use 
of other approximations (l)-(4), this replacement makes Eqs ( |2.13| )-( |2.15| ) a closed set of 
equations for a plasma including a liquid metal as a special case. The electron-ion and the 
ion-ion RDF's of liquid Rb at 313 K are plotted in Fig. ^ along with the effective ion-ion 
interactions calculated using both the QHNC and Geldart-Vosko G{Q); the resulting two 
effective potentials differ with each other, but yield almost the same ion-ion RDF's as shown 
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by the full curve and open circles. In Fig. |] the electron-ion RDF obtained from the QHNC 
method has an inner-core structure which is caused by the orthogonality of the free-electron 
to the bound-electron wave functions. On the other hand, it should be noticed that the usual 
liquid-metal theory based on the Ashcroft pseudopotential ||29| yields an electron-ion RDF, 
which has no inner-core structure (shown by full circles). The Ashcroft pseudopotential is 
constructed by neglect of this inner-core structure; this cut-off of the inner-core structure 
brings about a simple treatment of liquid metals in the standard liquid-metal theory. 

III. NUMERICAL CALCULATION APPLIED TO RUBIDIUM PLASMA 

Already, we have calculated the electronic and ionic structures of liquid rubidium in a 



wide range of temperatures and densities: compressed states ||3^ and expanded states |^ 
To show the applicability of the QHNC method to a plasma state, therefore here we take, 
as an example, liquid rubidium changing its temperature from 313 K to 3.5x10^ K (30 eV) 
at fixed ion-density rg = 5.388. Here, rl denotes the Wigner-Seitz radius a = rlaB in units of 
the Bohr radius qb- 

In the application of the QHNC equation to a plasma, we must calculate the electron- 
density distribution n'^{r\Ucs) under the external potential Ucs{r) to get gei{r) and ne(r|e) 
by solving the wave equation generally. However, it is difficult to determine this density 
distribution from the wave equation at high temperatures, where it becomes nearly the clas- 
sical Boltzmann factor: Tig(r|f/eff) = Sj/(ej)|^j(r)p =^ nQexp[—j3Ucs]. In the calculation 
of the free-electron density distribution, n^f{r\Ues) = 2/ /(ep)|-?/'p(r)p(ip/(27r^)^, at finite 
temperature, the electron kinetic energy ep is not limited within the Fermi energy Ep, as is 
the case at zero temperature. In addition, since the electron-density distribution n'^{r\Ueff) 
begins to approach the Boltzmann factor from the large distance as the temperature in- 
creases, we must calculate the wave functions with large angular momentum / to correctly 
obtain the classical electron-density distribution in the large distance at high temperature. 
We can circumvent this difficulty by using the Thomas-Fermi (TF) approximation to the 
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electron-density distribution for r > re] the distance Vc can be chosen according to the tem- 
perature so that n'^{r\Ucs) calculated from the wave equation becomes almost equal to the 
TF result for r > Vc- This situation can be seen in the following calculation. 

The electron-density distribution around a fixed electron is calculated for a partially 
degenerate electron plasma at the density 2.51 x 10^^/cm^ , that is, Ts = 4 in terms of 
rsflB = (3/47rnQ)^/^ defined for the electron density Uq; the results calculated for tempera- 
tures from 0.05 eV to 30 eV are shown in Fig. ^ The electron degeneracy is denoted by 
6 = k-QT/ E-pi i.e., the temperature over the Fermi energy. For high degeneracy (0.05 eV), 
the TF approximation (denoted by full circles) gives quite a different density distribution 
from the one calculated by the wave equation. When the temperature is increased to 10 eV, 
the TF result becomes almost same to that obtained by the wave equation except near the 
origin. When the temperature approaches 30 eV, the electron-electron correlation reduces 
to the classical one; the classical HNC equation for the one-component plasma (OCP) with 
the electron plasma parameter Fg = 0.277 (30 eV) provides an indistinguishable result from 
the TF calculation as is shown in Fig. |^ by the dashed curve. Here, the electron plasma 
parameter is defined by Fg = j3e^ /r^a-Q. In the calculation of ne(r|e), we can obtain at the 
same time the electron-electron DCF, which determines the LFC G{Q) from the relation 
j3Cee{Q) = —Vcc{Q)[i — GiQ)]'i the calculated G{Q) at r^ = 4 is shown in Fig. || for each 
temperature corresponding to Fig. ^. This figure indicates that the LFC at high temper- 
ature can be approximated by that calculated from the classical HNC for the OCP, as is 
shown by the case of 30 eV {6 = 9.62), where the LFC of the classical OCP is denoted by 
the dashed curve. In this way, by means of the QHNC equation for an electron gas we can 
obtain the electron-electron DCF, which determines the plasma properties in terms of the 
QHNC G{Q). 

Using now the QHNC-LFC G{Q) instead of a Geldart-Vosko ^ type G{Q), we apply 
Eqs. ( p.l3| )-( p.l4| ) to rubidium at the fixed density of the normal liquid metal; the temper- 
ature has been varied from (313 K) to 30 eV in order to investigate how a liquid-metal 
becomes a plasma. For the purpose to obtain initial data of the ion-core structure and 
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the electron-ion correlation for the fully self-consistent QHNC method, we take the jellium- 
vacancy model as a first step. In this model the following two approximations are introduced 
in the expression for the electron-ion interaction ( |2.11| ), that is, 



PU:!{r) = (3v,,{r) - Y. j C,,{\v - v'\)n[[gi,{r') - l]dv' . (3.1) 

(1) The ion-ion RDF involved above is approximated by the step function: gii{r) =9{r-a), and 

(2) the electron-ion DCF, by a pure Coulomb force: Cei(r) = — /?t>gj(r). Then, the problem 
to determine the electron-ion RDF becomes identical with the problem to determine the 
electron density distribution around a fixed nucleus at the center of the spherical-vacancy 
in the jellium background. This model is essentially same to the INFERNO model 0J^ (or 
the ion-sphere model) introduced by Liberman and also to the neutral-pseudoatom model 
|[T5| proposed by Perrot. By solving these integral equations, we can obtain the electron-ion 
DCF, and the ion-ion effective interaction ( p.l5|) with the combined use of the electron- 



electron DCF determined by the one-component QHNC equation ( p.l7| ). As a consequence, 
we can obtain initial input data for the fully self-consistent QHNC equations ( p^.l3D -( p.l5D . 



After that, this set of integral equations in conjunction with Eqs. ( f^.lT] ) and ( ^.181 ) is solved 



iteratively by varying the bound-electron number Z^ = j^ pTo{r)dY until the self-consistent 
condition (|2.4|) , n°{r)=p\,{r), is fulfilled. 

In the iteration to solve the QHNC equation, we must evaluate also the free-electron 
density distribution around the nucleus at arbitrary temperature. It has been the standard 
method [^ , p^ ,p|,^,[7|,p!8| for the calculation of this free-electron density distribution to use 



the TF approximation in such a way: 

nlf{r\U)^2f [exp(-/3(pV2m+[/(r)-/i°)) + l]-i-^ = nf(r|f/). (3.2) 

In this expression, the integration of p is limited within the domain p^/2r7i-|-f/(r) > to define 
the free-electron part of the usual TF formula, n^^{r\U), which involves both the bound- 
and free-electron density distributions. In Fig. ^, the electron-ion RDF's at temperatures 
of 3.5 X 10^ K and 3.5 x 10^ K calculated by the QHNC method are shown in comparison 
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with the results from the TF approximation ( p.2|) . It should be noticed that the TF formula 
for the free-electron density distribution ( |3.2| ) is not a good approximation except for large 
distances even at a high temperature of 3.5x10^ K [6=5.21), while there the total (free and 
bound) density distribution can be fairly well described by the TF formula. Noting this fact, 
we can evaluate the free-electron density distribution function ng'^(r |f/) for large distances 
r > Tc by the value of n^^ {r\U) —n'^'' {r\U) using the TF density distribution n'^^(r |f/) and 
the bound-electron density distribution n^^{r\U) obtained by the wave equation; in the case 
3.5x10^ K, the cut-distance Vc can be chosen to be 0.9a as shown by the arrow in Fig. |^. 
In this way, we can determine the the electron- ion RDF i.e., n^^{r\U), by solving the wave 
equation only for small angular momentum /; it is enough to take the maximum angular 
momentum /max = 15 even for the high temperature region. 

The electron-ion RDF's at temperatures of 1 eV and 3 eV are plotted in Fig. ^ along with 
the ion-ion RDF's and the effective ion-ion interactions. As is shown in Fig. |^, the electron- 
ion RDF at a temperature of 313 K has distinct inner-core and outer-core parts. Even at a 
temperature of 1 eV as shown in Fig. ^, this clear distinction remains as characteristic for a 
liquid metal; the ionization is practically unity and the ion-ion effective interaction is almost 
the same structure as that of liquid metal at the normal condition, although the ion-ion 
correlation becomes weak because of the small plasma parameter T = Z^ Fg. Therefore, 
we can consider that rubidium remains as a liquid-metal state even at 1.16x10^ K (1 eV). 
Figure ^ shows that at a temperature of 3 eV the distinction between inner- and outer-core 
parts near 0.4a disappears and that the ionization, now 1.21, has become significant. Because 
of the disappearance of this distinction between the the inner- and outer-core structures, it 
is difficult to construct a pseudopotenital in a plasma state. This makes a contrast with a 
liquid-metal state, where a pseudopotential such as the Ashcroft potential can be used to 
set up an effective interaction between ions in a liquid metal. 

The temperature variation of the electron-ion RDF is shown in Fig |^ for a range from 
eV to 30 eV. The electron-ion correlation becomes stronger for temperatures up to 10 eV, and 
turns to become weaker from 10 eV to 30 eV; the distinction between the inner- and outer- 
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core parts near the point 0.4a is gradually disappearing with increasing temperature. On the 
other hand, the ion-ion RDF's are shown in Fig. ^for increasing temperatures from eV (313 
K) to 30 eV. Also, the effective ion-ion interactions generating the RDF's in Fig. ^are plotted 
in Fig. 1^, where the open circles denote the screened Debye potential exp(— r/De)Zi^e^/r at 
a temperature of 30 eV (re=0.277) with De = (47re^/9nQ)^/^. We can see that the ion-ion 
effective potential approach the screened Debye potential as the electron plasma parameter 
Fe becomes small with increasing temperature. The HNC equation for a one-component fluid 
with the screened Debye potential at a temperature of 30 eV provides the RDF gii'(r) in 
fair agreement with the result from the QHNC method for the electron-nucleus mixture; for 
example, the structure factor at zero wavenumber becomes S'n(0) = 0.18 from the screened 
Debye potential, which should be compared with the QHNC result, S'n(O) = 0.16. This 
fact suggests that in the high temperature region where S'n(O) becomes large, the ion- 
sphere model (the jellium- vacancy model) can be improved by using the approximation 
gn{r) ^ gf^* (r) instead of the step function. 

Figure |TD| shows the temperature variation of the outer-bound levels of an average ion 
in Rb plasma at the fixed ion-density of normal liquid metal. The 4s- and 4p-bound levels 
are plotted there corresponding to free atom, 0, 1, 3, 5, 10, 22, and 30 eV, respectively. As 
the temperature increases, the bound levels become deeper due to the decrease of bound- 
electron number, which makes the screening effect weak. At temperatures of 22 and 30 eV, 
new bound levels, 5s and 4d, appear. The occupation number /(cj) at the level e^ is written 
at each level line in Fig. |10|. The ionization variation Zi is shown in the top of Fig. |1^ as 



the temperature is increased. 



IV. THE SAHA EQUATION IN THE DENSITY-FUNCTIONAL THEORY 

The DF theory provides the exact electron-density distribution ne{r\U) in the nonuniform 
electron system caused by an external potential U{r); however, it is important to notice 

16 



that this exact density distribution ne{r\U) at finite temperature is only an average density 
distribution. Consider a nucleus with the atomic number Za fixed in an electron gas. The 
effective external potential Ues{r) defined by the DF theory gives an exact expression for the 
electron density distribution in terms of a density distribution of noninteracting electrons in 
the form: 

ne(r|N) = n°(r|[/eff) = n^(r|N) + n^ (r|N) , (4.1) 

which defines an average ion with the Zb bound-electrons: Zb = J n]^{r\'N)dr. In the realistic 
system, the ion should have some integer-number of the bound-electrons with fiuctuations in 
time; the number Z^ given by the DF theory is only an average value of this bound-electron 
number over time. The similar situation is found in the ion structure in a nucleus-electron 
mixture. Here, we investigate the charge population of differently ionized species in a plasma 
on the basis of the DF theory. 

As is discussed in §. |T|, the average bound-electron number Z^ is defined by Eq. ( p.5| ) 
in the nucleus-electron model based on the DF theory, and the chemical potential /ig is 
determined by Eq. ( |2.7] ). This average bound-electron number Z^ in an ion in a plasma can 
be represented by 

dX ^exp[/?(ei-/xl^))] + 1 

if we introduce here the grand partition function SB = nf£i[l+'^6xp(— /9ej)]^* for the ion with 
bound-electrons, which have M bound-levels e^ with the degeneracy gi, and A = exp(/5/ig). 
Alternatively, the grand partition function Sb can be expanded in a polynomial of A: 

M G 

Sb = n[l + Aexp(-/36.)]^« = Y. ^^'^Zq (4.3) 

i=l (9=0 

with G* = I]^i Qi- In this expression, the canonical partition function Zq of the ion with the 
Q bound electrons is defined by 

ZQ^^fi(E,«)exp[-/3E«] (4.4) 
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M 

= E nU^n.exp(-/3e,n,)] (4.6) 

with the total energy -E^ = Z),tii ^s^^s = -E'^ \ {ns=0 or 1) for the Q bound-electrons. Here, 
il{Ei) represents the number of basic states with the energy Ei for the Q bound-electrons 
in the ion. Furthermore, the grand partition function Hb is written as Sb = J2f=i Uq using 
the function Uq defined by, 

M 

UQ^e^^''^ZQ= Y: Y[{ssCn.exp[/3ifi'i-es)ns]}. (4.7) 



From Eqs. ( [4.3|) and ( |4.7] ), we obtain the average bound-electron number Zb in another 
form: 



ZB = A-^lnSB= X^g^=<g>, (4.8) 

dA Q^o ^B 



which means that the probability for the bound-electron number of the ion to be Q is given 
by Uq/Eb. 

From another point of view, let us count the number of ions with the Q bound-electrons 
in a plasma, which is denoted by Nq; this satisfies the relation Y,q=o Nq = Ni since the total 
ion number in the system is A^i. In terms of Nq, the average bound-electron number Zb is 
determined by another way: 



G AT 

EQ^ = ^B. (4.9) 

Q=0 



Ni 



The above two expressions of Eqs. ([4.9| ) and ( |4.8|) for the probability that the ion in a plasma 
has Q bound-electrons give rise to the relation: Nq/Ni = Uq/Eb, that is, 

Nq Ni , , 

Uq ^b 

Since the righ-side of this equation is independent of Q, we obtain the expression: Nq/Uq = 
Nq^i/Uq-i, which can be rewritten in the form: 



Nq/Nq^, = Uq/Uq^i = expMZQ/ZQ_i . (4.11) 

If we introduce the canonical partition function Z'q = explPE^jZg using excitation energies 
E^—Eq measured from the ground state Eq of the ion with the Q bound-electrons, Eq. (|4.11|) 
is rewritten as 

^" exp[/3(/i° + /Q)]^ (4.12) 



with the ionization energy Iq = Eq — -Eq i^ the ground state. This is the Saha equation, 
which is applicable to a plasma where the electrons may be degenerate at any degree, and the 
ions and the electrons in a plasma may interact strongly with each other at high densities. 
When the temperature of a plasma becomes so high that the electrons behave as classical 
particles, the electron chemical potential is determined by the classical relation: 



R \3 '" ^f^ 



nlKll = ^ = exp[/3/x:;] (4.13) 

with the canonical partition function Z^, = 2(27rm//i^/3)'^/^ of a noninteracting electron gas 
and the thermal wavelength Ag. As a result from the above equation, Eq. (|4.12| ) is reduced 



to the usual expression for the Saha equation |34 



-^ = exp(/3/g)^^ (4.14) 

to determine the ion-density nq = Nq/V in the volume V. At this point, note that the 
fundamental relation ( ^.10| ) to derive the Saha equation is nothing but an ansatz introduced 



by Bar-Shalom et al. [B^ 



It should be recognized here that the solution of the Saha equation, ( ^4.12| ) or ( [4. 141 ) 



is obtained by determining the partition function Uq of a plasma with aids of the relation 
( [4.101 ), that is nq = UqUq/Eb- For this calculation, we can use the following recursive 
formulas I^S],^ for the partition function Uq of the ion with the Q bound-electrons: 



f/o = 1 , (4.15) 

Q 

UQ = J2XnUQ.n/Q, (4.16) 



n=l 
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where Xn = - E.^i ^^(-X^)" with X, = exp[-P{e,-fi% 

An apphed example of our formula to evaluate the ion population P{Q) is shown in the 
case of Rb plasmas. In the nucleus-electron model based on the DF theory, the bound-levels 
e, of the ion in a plasma is determined by solving the wave-equation for the self-consistent 
potential given by U^i{r) ( |3 . 1| ) , and the chemical potential /Xg of electrons is evaluated by 



the condition ( |2.7| ); the temperature variation of the bound-levels of Rb plasma was shown 



in Fig. [ly for a range from eV to 30 eV. Using these values, we can obtain Uq from the 



recursive relations, ( |4.15| ) and ( |4.16| ). In this way, the charge population P{Q) = Uq/Eb is 



evaluated for Rb plasma with the atomic number Za=37 varying the temperature from 3 
eV to 30 eV at fixed ion-density 1.03xl0^^/cm^ [r] = 5.388), as was studies in §. |TTl|for the 
evaluation of the electron-ion and ion-ion RDF's. Figure |ll| displays the charge population 
P{Q) in Rb plasmas for this temperature variation; the bound-electron number Z^ from 
Eq. ( [4.9|) is coincident with the values obtained previously by Eq. (|2.5| ) for each temperature, 
as a matter of course. At a sufficiently low temperature such as 1 eV, the charge population 
reduces to P{Q) = Sq^3q for Rb plasmas at liquid-metal density. 

V. CONCLUSIVE DISCUSSION 

We have demonstrated that the QHNC method, which has been successfully applied to 
many kinds of liquid metals, can be extended to treat a partially ionized plasma, taking 
rubidium as an illustrating example. In the application of the QHNC method to a plasma, 
it is necessary to use the LFC G{Q) at arbitrary temperature, which is determined by the 
one-component QHNC equation for an electron gas in the jellium; the QHNC G{Q) reduces 
to the LFC obtained by the classical HNC for the OCP at high temperatures as is shown 
in Fig. ^. In the numerical calculation in the QHNC method for a plasma, it is very time- 
consuming to evaluate free-electron density distributions at high temperatures; this problem 
can be easily circumvented by combined use of the TF approximation as discussed in §. p| . 
However, the simple TF approximation to the free-electron density distribution is shown to 
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give only a rough estimation in Fig. |^; this may not be apphed to calculate the accurate 
atomic structure in a plasma, although there are many examples as mentioned before. 

In a liquid metal, the electron-ion RDF has a clear inner-core structure distinct from 
outer-core structure (cf. Fig. ^; this distinction enables us to construct a pseudopotential 
in a liquid liquid metal. In the case of Rb, it remains as a liquid metal even at a temperature 
of 1.16 X 10^ K where the inner- and outer-core distinction is clearly seen with Z\ = \ as is 
displayed in Fig. |; this distinction disappears at 3 eV, where rubidium becomes a plasma 
with a significant ionization Z\ = \.2\. 

It is important to remember that the ion-sphere (IS) model is not an appropriate ap- 
proximation to treat a high-temperature plasma with a small plasma-parameter F. In the IS 
model, the ion-ion RDF is approximated by the step function with the Wigner-Seitz radius 
a. Therefore, this approximation is only valid in the strongly correlated region, where the 
structure factor S'ii(O) ~0 at zero wavenumber, because of the relation: 

- n\ Jlguir) - l]dr = 1 - 5n(0) ^ nf,47raV3 = 1 . (5.1) 

This condition is very important to keep charge neutrality around the ion P7| : 

Zi = -Zml Jlguir) - l]dr + n^ j [g,i{r) - l]dv . (5.2) 

When the ion-ion RDF becomes weaker as the temperature increases, the structure factor 
5*11(0) grows large; there, the IS model is not properly applicable. This situation is exem- 
plified in our calculation of the ionization Zi. The IS (jellium- vacancy) model provides the 
ionization Zj: 1.27, 2.05, 3.80, 5.28 and 6.55 for temperatures: 3, 5, 10, 22, and 30 eV, 
respectively, while the corresponding values from the fully self-consistent QHNC calculation 
are 1.21, 1.96, 3.71, 5.08 and 6.08, respectively. We can see that significant differences are 
manifested between the IS and QHNC results, as the temperature increases. The reason for 
these differences can be ascribed to the fact that the IS model produces a screening effect 
from the ions contained in Eq. ( |5.2| ) in the approximated form Zi[l — 5*11(0)] ~ Zi, which 
becomes too strong for a plasma at high temperatures where 5*11(0) becomes large. As a 
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consequence, the neutral-atom model [^] based on the IS model is, also, not appropriate 
to construct an effective ion-ion interaction at a high-temperature plasma. As discussed in 
§.|rn|, the IS model can be improved by using the ion- ion RDF from the HNC equation for 
a screened Debye potential instead of the step function, in dealing with a high-temperature 
plasma. 

In the QHNC method, the inner-structure (the atomic structure of the ion) is determined 
in the consistent way with the outer structures (the electron-ion and ion-ion RDF's and 
the average ionization Zi). Therefore, we can expect this method to provide an accurate 
procedure to deal with the atomic structure in a high-density plasma; the bound-levels in 
an ion can be calculated by taking account of the density and temperature effects as is 



shown in Fig. 10. In addition, it should be remarked that the DF theory leads to the 



Saha equation as discussed in §. |I^, and the QHNC method based on the DF theory can 
provide a procedure to solve the Saha equation with ease by using the recursive formula. 
As an applied example of this formula, the charge population P{Q) is calculated from the 
QHNC result for a Rb plasma, as is displayed in Fig. ^. Moreover, we can expect that the 
QHNC method can solve various kinds of problems associated with the atomic structure in 
a plasma by taking account of the plasma effects. For example, with the combined use of 
Slater's transition-state method [^, we have already calculated the shift-variation of the 



K-edge |p9| in an aluminum plasma along the shock Hugoniot in good agreement with the 



experiment performed by DaSilva et al. |40 



The QHNC method can provide an accurate description of the metallic system for a wide 
range of densities and temperatures from the liquid-metal to the plasma state in a unified 
manner, as is ascertained from many experiments on liquid metals. This fact indicates that 
the QHNC method can be used to calculate transport properties and equation of state in a 
wide region from the liquid-metallic to the plasma state, where there has been no systematic 
applicable theory up to the present. 

With decreasing temperature or increasing pressure of a plasma, some bound state of 
each ion in a liquid metal or plasma begins to disappear into the continuum; it becomes 
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a narrow-resonant state and disappears gradually as a wide resonance in the continuum. 
In our calculation of plasma states we do not take account of resonant states. In practice, 
the resonant-state contribution in a plasma is not as significant as in the case of a liquid 



metal such as a transition metal. A precise definition of a resonant state |^ is given by 
the pole Ene of the S-matrix Si{E) concerning the wave equation for an electron under the 
effective potential ( p.2|) with Se{E) = exp{2i6i{E)) for phase shifts Si{E). In a strict way, it 
is required ^ that the "bound" -electron number Zb in an ion should include a contribution 
of the physical resonant states ( \QmEne\ -C ^eEnf, ) in addition to the bound electrons with 
ej < in such a way that 

Zb=Y. f{e,) + Y. 2(2£ + me[F{K,)] . (5.3) 

ei<0 ?i£gphys.res. 



In the above, F{Ene) is the function introduced by More |0 with the definition: 



to represent the thermal occupation probability of a resonant state Eni. Also, the chemical 
potential /ig in consideration of the resonant states should be determined by 

nlJ exp[/5(p2/2m - /i^)] + 1 {2TTn)^ " ^ ' ' 

However, the determination of the charge occupation P{Q) taking account of the resonant- 
state contribution is a problem which remains to be investigated. 

Molecular-dynamics (MD) simulation is necessary to study complex systems which are 
inhomogeneous or time dependent and so on. In the MD simulation of a dense plasma, the 
Coulomb interactions among close and distant particles must be calculated precisely and 
efficiently; the Particle-Particle Particle-Mesh (PPPM) method [^,0] should be used in the 



simulation code to treat many particles. In the particle-particle method the Coulomb forces 
among close particles are directly summed up and in the particle-mesh method the forces 
on a particle are interpolated from electric fields at mesh points. In the SCOPE (Strongly 
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COupled Plasma ParticlE) code [^^-|46[] based on the PPPM method, the Deutsch potential 



P7| is adopted to imitate quantum effects. With use of the code, bremsstrahlung emission 



45| , transport coefficients and the Lyapunov exponents [^ were obtained in dense plasmas. 
However, the applicability of the Deutsch potential is limited to a hydrogen plasma or to a 
fully ionized plasma at most, since the ion structure is not considered in the derivation of the 
Deutsch potential. In order to perform a classical MD simulation (SCOPE) on a partially 
ionized plasma, we must introduce effective classical potentials v^j (r) applicable to partially 
ionized ions in a plasma as a classical mixture of ions and electrons; the QHNC method 
can produce these effective potentials as follows. The quantum effects of electron-electron 
interaction can be taken into account by defining an effective classical pair-potential v^^{r) 
between electrons in such a way that the HNC equation for nl{r\e) = nQg^^{r) in classical 
fluids with fee(^) provides the same electron-density distribution n^^^'~' {r\e) determined by 



the one-component QHNC equation (|2.17|) ; this condition is written in the following integral 
equation for v^^{r), 

nl{r\e) = n^exp[-/3<,(r)+7^(r)] = nf^^^{r\e) (5.6) 

with 7'^(r) = / Cpg(|r— r'|)[r2g(r'|e)— r2Q](ir' . In a similar way, an electron-ion classical potential 
fgj(r) is determined by the condition that the classical electron-ion RDF g^iir) should be 
identical with the QHNC result: 

g:,{r) ^ exp[-/?<i(r) + r:,{r)] = <?,™(^) (5-7) 

with 

^elW = /Qe(|r-r1)n^o[^:i(^')-l]rfr' + |Q(|r-r'|)ni[^if^^(r')-l]rfr. (5.8) 

With use of the effective potentials determined above, the SCOPE code can be applied 
to investigate dynamical problems in a partially ionized plasma as a classical ion-electron 
mixture. 

We have shown that the QHNC method is extended to treat a partially ionized plasma 
in a wide range of densities and temperatures, and provides the average ionization, the 
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electron-ion and ion-ion RDF's, the atomic structure of the ions and the charge population 
of differently ionized species in a self-consistent manner from the atomic number as the 
only input data. Therefore, this method produces the fundamental quantities necessary to 
calculate the plasma properties, and offers a procedure to treat the spectroscopic problem 
in a plasma. It should be kept in mind that the QHNC method can provide a precise 
description of "simple" plasma where the bound states are clearly distinguished from the 
continuum state; to take into account the resonant states in a plasma, some improvement is 
necessary as was discussed in the previous work |]l|. 
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FIGURES 
FIG. 1. Ion-ion structure factor Sii{Q) for liquid Rb at a temperature of 313 K; the QHNC 

method yields a structure factor (full curve) in excellent agreement with experiments (open ||2^ 

and full circles 



FIG. 2. The electron-ion and ion-ion RDF's with the effective interactions in liquid Rb at 313 
K. The solid curves denote the RDF's, the total electron-density distribution and the effective 
ion-ion interaction calculated by using the QHNC G{Q). The dashed curve designates the effective 
ion- ion potential based on the Geldart-Vosko G{Q), which yields the ion-ion RDF plotted by the 
open circles; •, the electron-ion RDF derived by using an Ashcroft potential. 

FIG. 3. Electron-electron correlations in an electron gas at the density of rg = 4 for temperatures 
ranging from 0.05 to 30 eV with 6 = k-QT/E-p indicating the electron degeneracy. The solid curves are 
calculated by the one-component QHNC equation, while the full circles denote the result from the 
TF approximation for each temperature. The classical HNC equation for the OCP corresponding 
to a temperature of 30 eV provides the electron-electron correlation plotted by the dashed curve 
("OCP"), which is indistinguishable from the TF result at a temperature of 30 eV. 

FIG. 4. The local-field corrections (LFC) at temperatures from to 30 eV determined by the 
one-component QHNC equation. The LFC calculated by the classical HNC for the OCP with 
Fe = 0.227 is displayed by the dashed curve, which indicates that this can be used to approximate 
the QHNC LFC at this temperature. 
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FIG. 5. The electron-ion RDF's calculated from the wave equation and the TF approximation 
at temperatures of 3.5 x 10^ K (3 eV) and 3.5 x 10^ K (30 eV). The solid curves denote the results 
from the wave equation and the dashed curves, from the TF approximation. The TF approximation 
can not give a good description of the free-electron density distribution ng(r[I) =nQgci(r) in the 
core-region even at a high temperature (30 eV), where the total (bound and free) electron-density 
distribution is fairly well described by the TF approximation. The electron-ion RDF can be 
approximated by the TF formula for larger distances than the point (0.9) denoted by the arrow. 

FIG. 6. The electron-ion and ion-ion RDF's together with the effective ion-ion interaction at 
temperatures of leV and 3 eV. The electron-ion RDF at 1 eV has an inner-core structure similar 
to the RDF of a liquid state at the normal condition (cf. Fig. |^, and this inner-core structure near 
r/a = 0.4 disappears at 3 eV with a significant ionization Zi = 1.21. 

FIG. 7. The temperature dependence of the electron-ion RDF for a range from to 30 eV. The 
inner-core structure near 0.12 reflects the variation of the bound-electron wave functions in an ion 
due to the orthogonality between the free- and bound-electron wave functions. 

FIG. 8. The temperature dependence of the ion-ion RDF for a range from eV (313 K) to 30 
eV. The ion-ion RDF becomes small at 1 eV; nevertheless rubidium remains as a liquid state with 
Zi = l. 

FIG. 9. The temperature dependence of the effective ion- ion interaction for a range from eV 
(313 K) to 30 eV. The effective potential at a temperature of 30 eV approaches the screened Debye 
potential denoted by the open circles. 
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FIG. 10. The temperature variation of outer bound levels (4s, 4p, 4d and 5s) in the Rb ion in 
a plasma for a range from eV to 30 eV at a fixed density of rl = 5.388. The bound levels become 
shallow when free atoms are compressed to be a liquid state (0 eV), and turn to become deeper as 
the temperature is increased with the fixed ion-density. Numbers attached to bound levels denote 
the occupation numbers /(cj), and the ionization Zi for each temperature is written in the top of 
this figure. 

FIG. 11. The dependence of the charge population P{Q) on the temperature varying from 3 
eV to 30 eV for Rb plasma (2'a = 37) at a fixed ion-density of 1.03 x 10^2 /cm^ {rl = 5.388). This 
charge population provides the average bound-electron Zb of ion, as denoted in this figure for each 
temperature along with the ionization Zi = Za — .^b- 
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Figures Figures 




Q (A-^) 

FIG. 1. Ion-ion structure factor Sii{Q) for liquid Rb at a 
temperature of 313 K; the QHNC method yields a structure 
factor (full curve) in excellent agreement with experiments 
(open [26] and full circles [27]). 
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FIG. 2. The electron-ion and ion-ion RDF's with the ef- 
fective interactions in liquid Rb at 313 K. The solid curves 
denote the RDF's, the total electron-density distribution and 
the effective ion-ion interaction calculated by using the QHNC 
G{Q). The dashed curve designates the effective ion-ion po- 
tential based on the Geldart-Vosko G{Q), which yields the 
ion-ion RDF plotted by the open circles; •, the electron-ion 
RDF derived by using an Ashcroft potential. 
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FIG. 3. Flectron-electron correlations in an electron gas 
at the density of fg = 4 for temperatures ranging from 0.05 
to 30 eV with 9 = /cbT/Bf indicating the electron degener- 
acy. The solid curves are calculated by the one-component 
QHNC equation, while the full circles denote the result from 
the TF approximation for each temperature. The classical 
HNC equation for the OCP corresponding to a temperature 
of 30 eV provides the electron-electron correlation plotted by 
the dashed curve ("OCP"), which is indistinguishable from 
the TF result at a temperature of 30 eV. 
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FIG. 4. The local-field corrections (LFC) at temperatures 
from to 30 eV determined by the one-component QHNC 
equation. The LFC calculated by the classical HNC for the 
OCP with Fe = 0.227 is displayed by the dashed curve, which 
indicates that this can be used to approximate the QHNC 
LFC at this temperature. 
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FIG. 5. The electron-ion RDF's calculated from the wave 
equation and the TF approximation at temperatures of 
3.5 X 10'' K (3 eV) and 3.5 x 10^ K (30 eV). The soUd curves 
denote the results from the wave equation and the dashed 
curves, from the TF approximation. The TF approximation 
can not give a good description of the free-electron density 
distribution ne(r|l) = nogei(r) in the core-region even at a 
high temperature (30 eV), where the total (bound and free) 
electron-density distribution is fairly well described by the TF 
approximation. The electron-ion RDF can be approximated 
by the TF formula for larger distances than the point (0.9) 
denoted by the arrow. 
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FIG. 6. The electron-ion and ion-ion RDF's together with 
the effective ion-ion interaction at temperatures of leV and 3 
eV. The electron-ion RDF at 1 eV has an inner-core structure 
similar to the RDF of a liquid state at the normal condition 
(cf. Fig. 2), and this inner-core structure near r/a = 0.4 dis- 
appears at 3 eV with a significant ionization ^1 = 1.21. 
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FIG. 7. The temperature dependence of the electron-ion 
RDF for a range from to 30 eV. The inner-core structure 
near 0.12 reflects the variation of the bound-electron wave 
functions in an ion due to the orthogonality between the free- 
and bound-electron wave functions. 
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FIG. 9. The temperature dependence of the effective 
ion-ion interaction for a range from eV (313 K) to 30 eV. 
The effective potential at a temperature of 30 eV approaches 
the screened Debye potential denoted by the open circles. 



FIG. 8. The temperature dependence of the ion-ion RDF 
for a range from eV (313 K) to 30 eV. The ion-ion RDF 
becomes small at 1 eV; nevertheless rubidium remains as a 
liquid state with Zi = l. 
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FIG. 10. The temperature variation of outer bound levels 
(4s, 4p, 4d and 5s) in the Rb ion in a plasma for a range 
from eV to 30 eV at the fixed density of fg = 5.388. The 
bound levels become shallow when free atoms are compressed 
to be a liquid state (0 eV), and turn to become deeper as the 
temperature is increased with the fixed density. Numbers at- 
tached to bound levels denote the occupation numbers f{^i), 
and the ionization Zi for each temperature is written in the 
top of this figure. 
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FIG. 11. The dependence of the charge population P{Q) 
on the temperature varying from 3 eV to 30 eV for Rb plasma 
(Za=37) at the fixed density of 1.03xl0^Vcm^ (r,^ = 5.388). 
This charge population provides the average bound-electron 
Zb of ion, as denoted in this figure for each temperature along 
with the ionization Zi = Za — Zb • 
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